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Abstract Solar arrays are structures which are connected to satellites; during launch, 
they are in a folded position and submitted to high vibrations. In order to save mass, 
the flexibility of the panels is not negligible and they may strike each other; this may 
damage the structure. To prevent this, rubber snubbers are mounted at well chosen 
points of the structure; a prestress is applied to the snubber; but it is quite difficult to 
check the amount of prestress and the snubber may act only on one side; they will 
be modeled as one sided springs (see figure|2]i. 

In this article, some analysis for responses (displacements) in both time and fre- 
quency domains for a clamped-clamped Euler-Bernoulli beam model with a spring 
are presented. This spring can be unilateral or bilateral fixed at a point. The mount- 
ing (beam +spring) is fixed on a rigid support which has a sinusoidal motion of 
constant frequency. 

The system is also studied in the frequency domain by sweeping frequencies be- 
tween two fixed values, in order to save the maximum of displacements correspond- 
ing to each frequency. Numerical results are compared with exact solutions in par- 
ticular cases which already exist in the literature. 

On the other hand, a numerical and theoretical investigation of nonlinear normal 
mode (NNM) can be a new method to describe nonlinear behaviors, this work is in 
progress. 
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Fig. 1 At left: solar arrays from folded to final position, at right: picture for the satellite AMC12 
from Thales Alenia Space company 



1 Introduction 

1.1 Previous Works 

In articles (TJ, JfQ, J5) et iflOl . proposed by Thales Alenia Space research team, the 
dynamic of a beam sytem with a nonlinear contact force, under a periodic excitation 
given as an imposed acceleration form is studied both numerically and experimen- 
tally. When sweeping frequencies in an interval which contains eigen frequencies 
of the beam, resonance phenomena appear as well as new frequencies caused by the 
unilateral contact. 

Finite element method in space domain is used, followed by numerical integration 
of the ordinary differential systems using specific software like 'STRDYN' of the 
finite element package DIANA. 

The frequency sweeping is done in different ways, one of these way is such as the 
frequency / changes as a function of the time t according to :f(t ) = fo2 st / 60 where 
s is the sweep rate in octaves/min and /o is the start frequency of the sweep. Results 
prove differences between sweep-down and sweep-up around eigen frequencies of 
the system where solutions are unstable. 

At each value of the time t , computation is done and the maximum of displacement 
and acceleration are saved. 

Finally, comparison is made in time and frequency domains, between linear and 
nonlinear cases. 
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Fig. 2 left: Simplified Mechanical Model, right: folded solar arrays with snubbers 



1.2 Present Work 

This work is part of the phd work of the author under the guidance of B. Rousselet 
with the support of Thales Alenia Space, France. 

Some analysis for responses (displacements) in both time and frequency domain 
for a clamped-clamped Euler-Bernoulli beam model with a linear spring are pre- 
sented, the spring can be unilateral or bilateral fixed at a point. The mounting (beam 
+spring) is fixed on a rigid support which has a sinusoidal motion of constant fre- 
quency. 

The system is also studied in the frequency domain by sweeping frequencies be- 
tween two fixed values, such as saving the maximum of displacements correspond- 
ing to each frequency. 

Numerical results are compared with exact solutions in particular cases which 
already exist in the literature. 

On the other hand, a numerical and theoretical investigation of nonlinear normal 
mode (NNM) can be a new method to distinguish linear from nonlinear cases. 



2 Simplified Mechanical Model 

The study of the total dynamic behavior of solar arrays in a folded position with 
snubbers are so complicated, that to simplify, a solar array is modeled by a clamped- 
clamped Bernoulli beam with one-sided linear spring. This system is fixed on a 
shaker which has a vibratory motion d(t) see figure (ff]). 

The motion of this beam system is modeled by the following PDEs with bound- 
ary conditions : u(Q,t) = u(L,t) = d(t) and d x u(Q,t) = d x u(L,t) = 0: 
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Bilateral spring: pSu(x,t) +EIw-"\x,t) = k,-(d(t) - u(x ,t))8 X0 (1) 
Unilateral spring: pSii(x,t)+EIu^' v \x,t) = k,-(d(t) - u(x Q ,t)) + 8 XQ (2) 

u + is the function denned by u + = L= 0.485 m beam length. 

k r = spring stiffness. d(t) = — " sin(27r/f) is the shaker motion. 

p = 2700 kg/m 3 beam density. S = 7,5. 10~ 4 m 2 cross sectional area. 

E = 7. 10 10 N/m 2 Young's modulus. 1, 4 1. 1 0~ 8 m 4 second moment of area. 

This data are taken from [4 a . The classical Hermite cubic finite element approxima- 
tion is used here to find an approximate solution for equations (Q~} and ©, we find 
then two ordinary differential systems in the form : 

Mq + Kq = k r (d(t) - q XQ )e? (3) 

Mq + Kq =k T {d{t)- q X(l )+e? (4) 

M et K are respectively the mass and the stiffness assembled matrices, q is the vec- 
tor of degree of freedom of the beam, qi ~ (iij,d x Ui), i = 1,2, ...,n. To each node are 
associated two degrees of freedom, the displacement and its derivative. 

To integrate numerically systems ([3]) and Q, we use the Scilab routines "ODE's" 
followed by the FFT (Fast Fourier Transformation) to find frequencies of solu- 
tions, there is no special treatment for "ODE" routine to deal with the local non- 
differentiable nonlinearity (d(t) — u(xo,t)) + . 



3 Numerical Results 



3.1 Highlights on the linearized system and the nonlinear effects 

There are two linear cases, the first one when there is no spring attached to the 
beam, and the second case is such as a linear spring is always attached to the beam 
( bilateral spring ). The system becomes nonlinear when the spring becomes one- 
sided when the prestress is not well tuned. The nonlinearity has a special form, it is 
locally not differentiable but it is lipchitz. 

Without spring, the motion has the following linear system of equation : 

Mq + Kq = (5) 

The eigen frequencies of the motion are easily calculated by computing the gener- 
alized eigen values of M and K, this calculus can be done using software like Scilab 
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j9]- We are interested by the first three eigen frequencies of the system, we note 
that the precision depends on the number of the finite elements used for the mod- 
eling of the system. We verify that ten finite elements give a good approximation 
for the first three eigen frequencies and their values are respectively :334. 21622Hz, 
921.48815Hz and 1807.7966Hz. As known, the eigen-frequencies of a Bernoulli 
beam could be calculated using this formula: 

/'■ = T&yjtimt' V* i = 1,2,3, ...» are given in @. m = 4.73, ^2 = 7.853 and 
H 3 = 10.996 give ft = 334.19889Hz, ft = 921.19996Hz and ft = 1806.1432Hz. 
The nonlinearity on the system modifies the motion by adding new frequencies, 
subharmonics and superharmonics, besides the eigen frequencies. When the system 
is under periodic excitation, new frequencies appear also, there are many combina- 
tions of the excitation frequency with the system frequencies. These new frequencies 
appears in FFT of the system and also in the sweep test. There are many ways to 
calculate these frequencies, harmonic balance method and nonlinear normal mode 
(MNN) and asymptotic expansions methods (|7|) . 



3.2 One Node Finite Element Model Without Periodic Excitation 

In this case, the beam is modeled by two finite elements without periodic excitation. 
Equations of displacement and its derivative are independent here because of the 
structure of the mass and stiffness matrices: 

/ 0.3647893 \ A 66 1090 \ 

M ~y o 0.0005500 J ~ \ 32560.825 )' 

To find eigen frequencies of the system in the linear case (system without contact), 
we juste have to calculate the generalized eigen values of M and K. 
The motion is divided into two phases: the first one when the beam touches the 
spring and another one when the beam does not touch the spring. The spring mass is 
negligible beside the beam mass, the equations of displacement of these two phases 
are respectively : 

M(\,\)u+{K(\,\)+k r )u = 

M(l,l)ii + K(l,l)u = 

The boundary conditions are :u(Q,t) = u(L,t) = and d x u(0,t) = d x u(L,t) = 0. 
The period of the solution will be the sum of the half period of the first phase and 
the second phase : 

K % 

T = — ^=^= H 

/ Kll(U) / A: n (l,l)+A> 
Vm u (i,i) V M„(l,l) 

the numerical value of the motion frequency j is 384. 74186Hz this is in agreement 
with the numerical calculus. In figure Q, The Fast Fourier Transformation shows 
frequencies of the system, the first peak correspond to the analytical value. The 
motion in the phase plane shows a periodic conservative solution. 
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Fig. 3 Left: FIT of the displacement of the second node with unilateral spring contact with initial 
conditions u(x,0) = 0.4, u(x, 0) = 0. Right: solution in the phase space. 

3.3 One Node Finite Element Model, Beam Under Periodic 
Excitation 

The beam is modeled by two finite elements under periodic prescribed displacement 
d(t). The spring is always fixed in the middle of the beam on a node, it can be uni- 
lateral or bilateral. 

Figure (0]i shows the displacement of the second node (in the middle) of the beam 
with bilateral and unilateral spring, the period of the motion changes with unilateral 
spring. Figure (O shows frequencies of the beam with bilateral and unilateral cases; 
In the bilateral case, there are two frequencies, the eigen frequency (430 Hz) and 
the frequency of the shaker (500 Hz); In the unilateral case, there are the system 
frequency, the shaker frequency and many other superharmonics and subharmonics 
frequencies due to the contact. 



3.4 Ten Finite Elements 

In this section, the beam is modeled with ten finite elements, this approximation is 
quite good if we are interested by the first two eigen frequencies of the system. 
The spring is always fixed in the middle of the beam on a node, and can be unilateral 
or bilateral, the whole system is fixed on a shaker which has periodic prescribed 
displacement d{t), we compare the linear with the nonlinear case for the same value 
of parameters. 

Figure (O shows the displacement in the time domain for bilateral and unilateral 
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case and their corresponding frequencies using FFT. In the linear case (bilateral 
spring), the first and the third eigen frequencies are shown, the shaker frequency 
also(500 Hz); The middle of the beam is node of the second mode. In the nonlinear 
case (unilateral spring), frequencies corresponding to the first and the third eigen 
frequencies of the linear case and also the shaker frequency are shown; Many other 
superharmonics, subharmonics and combination with the shaker frequency appear 
also due to the non linear contact. 



4 Frequency Sweep Excitation 

Frequency sweep excitation is usually used in the experiments to check the dynamic 
behavior of mechanical systems. As mentioned before, the beam is under periodic 
excitation given as a displacement form 

d(t) = — [2k f) 1 sm (^ 7i: ft)> mat means that the beam is under an effort (acceleration) 
d(t) = asm(2nft), a in m/s 2 is the amplitude of the acceleration and / in Hz is the 
frequency of excitation. 

We are just interested by a sweep-up test, for a initial given frequency fo and a 
fixed value of a. We compute the solution of differential systems, then we save the 
maximum of the acceleration and the displacement. In a second time, we add a fixed 
value to /o, the frequency step df. We compute again for fo +df and we save the 
same quantities for this iteration, the initial conditions are fixed on zero again. We 
continue our test to reach a fixed frequency f\ chosen such as to cover the first and 
the second eigen frequency of the system. 

Finally, we plot the maximum saved in each sweep test, then we compare linear 
with nonlinear cases, we study also the effect of the amplitude a and the frequency 
step on the system. 

We remark that obtained curves are similar to curves of FFT, it is an another 
method to find frequencies for mechanical systems. 



4.1 Two Finite Elements 

The beam is modeled by two finite elements, the spring is fixed in the middle, the 
Bernoulli beam is clamped in its both extremities, we just have a free node of two 
degrees of freedom, the first is the displacement, the second is the derivative of the 
displacement. The eigen frequency of the linearized system in this case is around 
339Hz, the sweep-up begin from 100Hz to 1000Hz, the frequency step df = 5Hz 
and the amplitude of excitation a = 5Qm/s 2 . 

In each iteration, the integration time is tf = 0.1s, the initial conditions are al- 
ways fixed at : q(x,0) = q(x,0) = 0,k r = l0 6 N/m. 
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In figure (0, the peak in the bilateral case corresponds to the eigen frequency of 
the system, its abscissa is 430//z. In the unilateral case, the peak corresponds to the 
eigen frequency too, its abscissa is 384f/z. Other peaks appear; they are due to the 
unilateral contact. These results conform well with the FFT in figure ©. 



4.2 Ten finite elements 

The beam is modeled by ten finite elements, the spring always in the middle, the 
other parameters are the same defined in section l4~T1 

In figure (O, the peak in the bilateral case corresponds to the first eigen frequency 
of the system, the second eigen frequency does not appear because the imposed 
displacement is not enough to excite it (d(t) = — n " a sin(2?r/f ) very small when 
f become so high). Peaks in the unilateral case show the system frequencies corre- 
sponding to the eigen frequency of the linear system and other frequencies due to the 
contact. These results conform with the FFT in figure © in the interval [100, 1000] 
Hz. 




Fig. 4 Displacement of the second node with bilateral and unilateral spring contact under peri- 
odic excitation of 500 Hz (Two finite element model) 
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Fig. 6 Displacement and their FFT of the sixth node with bilateral and unilateral spring contact 
under periodic excitation of 500 Hz (Ten finite element model) 
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Fig. 7 Maximum of absolute value of the displacement for both bilateral and unilateral spring, 

two finite element model, kr = l0 6 N/m, a — 50m/s 2 , tf = 0. Is (Sweep-up test) 
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Fig. 8 Maximum of absolute value of the displacement of all nodes(log scale) for both bilateral 
and unilateral spring, ten finite element model k r = 10 6 , tf = 0.1.S, a = 50m /s 1 (sweep-up test) 

5 Conclusion 

We have presented some preliminary numerical results to compare the vibrations of 
a beam equipped with a bilateral or a unilateral spring. Asymptotic expansions using 
some results of S. Junca ([7]) are in project and compared with numerical results in 
order to asses the quality of both approaches. In particular accurate computations of 
non linear normal modes (see [3]) will be considered. 
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